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ABSTRACT 

Predictions from a Galactic Structure and Kinematic model are compared to the 
absolute proper-motions of about 30,000 randomly selected stars with 9 < i?j < 19 
derived from the Southern Proper-Motion Program (SPM) toward the South Galactic 
Pole. The absolute nature of the SPM proper-motions allow us to measure not only the 
relative motion of the Sun with respect to the local disk, but also, and most importantly, 
the overall state of rotation of the local disk with respect to galaxies. The SPM data 
are best fit by models having a solar peculiar motion of +5 km s~^ in the V-component 
(pointing in the direction of Galactic rotation), a large LSR speed of 270 km s"^, and a 
disk velocity ellipsoid that points towards the Galactic center. We stress, however, that 
these results rest crucially on the assumptions of both axisymmetry and equilibrium 
dynamics. 

The absolute proper-motions in the U-component indicate a solar peculiar motion of 
11.0 ± 1.5 km s~^, with no need for a local expansion or contraction term. 

The implications of the large LSR speed are discussed in terms of gravitational mass 
of the Galaxy inferred from the most recent and accurate determination for the proper- 
motion of the LMC. We find that our derived value for the LSR is consistent both with 
the mass of the Galaxy inferred from the motion of the Clouds (3 — 4 x IO^^Mq to 
~ 50 kpc), as well as the timing argument, based on the binary motion of M31 and the 
Milky Way, and Leo I and the Milky Way (> 1.2 x 10^^^© to ~ 200 kpc). 

Subject headings: Astrometry: stellar dynamics - stars: kinematics - stars: fundamental 
parameters - Galaxy: fundamental parameters - Galaxy: kinematics and dynamics - 
Galaxy: structure 
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1. Introduction 

In this paper, we present the main rcsuhs from an 
analysis of the kinematic data obtained in the con- 
text of a massive absolute proper-motion survey, the 
Southern Proper-Motion Program (SPM hereafter). 
Our analysis is basc^d on the SPM Catalog as de- 
scribed by Platais ct al. (1998). This catalog provides 
positions, absolute proper-motions, and BV photom- 
etry for about 30,000 randomly selected stars, among 
other objects, near the South Galactic Pole. The sky 
coverage of the SPM Catalog is about 720 deg^ in 
the magnitude range 5 < V < 18.5. The accuracy 
of individual absolute proper- motions is 3-8 mas yr~^ 
depending on the star's magnitude. For a complete 
description of the catalogue structure, contents, plate 
measurement and other astrometric/photometric de- 
tails, the reader is referred to Girard et al. (1998, 
Paper I) and Platais et al. (1998, Paper II). A de- 
tailed analysis of the velocity distribution function, 
and star- and color-counts for stars from the disk, 
thick-disk, and halo derived from the SPM photom- 
etry and astrometry is addressed by Mendez et al. 
(1999, Paper III). 

The arrangement of this Letter is the following: In 
Section 2 we present an overview of the SPM pho- 
tometric and proper-motion data. In Section 3 the 
model used to compare the observed proper-motion 
distributions with those predicted is described, while 
Section 4 presents the model predictions as compared 
to the proper-motion data. Section 5 presents the 
implications of our results for the global mass of the 
Galaxy. Finally, Section 6 presents the main conclu- 
sions of the paper. 

2. The randomly selected sample 

The random SGP-SPM sample has been chosen 
such that, at a given Bj magnitude interval, a fixed 
number of stars are randomly extracted from all 
stars available in that magnitude interval. As a re- 
sult, one-magnitude intervals are chosen in the range 
9 < Bj < 19 for the subsequent analysis. At a given 
magnitude interval, the complete sample and the ran- 
domly selected sample differ only by a scale factor. 
Therefore, color and kinematic properties binned in 
the proper magnitude intervals should be very simi- 
lar to those derived from a complete sample, except 
for the increased uncertainties on the derived values 
because of the smaller sample - an effect that is fully 
taken into account in the analysis below. 



3. The starcounts cind kinematic model 

We model the starcounts concurrently with the 
kinematics by using the Galactic structure and kine- 
matic model presented by Mendez and van Altena 
(1996). The starcount model employed here has been 
tested under many different circumstances, and has 
proved to be able to predict starcounts that match 
the observed magnitude and color counts (in both 
shape and number) to better than 10%, and in many 
cases to better than 1%. The kinematic model pre- 
sented by Mendez and van Altena (1996) has also 
been shown to be able to reproduce the kinemat- 
ics of disk stars toward different Galactic fields. In 
Mendez et al. (1999), a full analysis of the SPM- 
SGP data in terms of star and color-counts, as well 
as the kinematics is presented. A basic assumption 
in the Mendez and van Altena (1996) model is that 
of the Gaussian shape of the velocity ellipsoid. While 
the velocity distribution observed for late-type stars 
(e.g., Dehnen 1998) is clearly skewed and cannot be 
properly modeled by a Gaussian, the error in the pre- 
dicted median proper motion made by using a shifted 
Gaussian is likely to be small. Indeed, Mendez et al. 
(1999), show that shifted- Gaussian velocity distribu- 
tion functions do provide a good fit to the full SPM 
proper-motion distributions as a fmiction of apparent 
magnitude, with perhaps a slight predicted excess of 
stars with positive proper-motions in the direction of 
Galactic rotation. 

Because of projection effects near the SGP, our 
model comparisons are most sensitive to the motions 
in the direction of Galactic rotation (the so-called V- 
component), and along the Galactic center-anticenter 
direction (the U-component), and not to the mo- 
tions perpendicular to the Galactic plane (the W- 
component). 

3.1. The Solar peculiar velocity and the mo- 
tion of the LSR 

There have been a number of determinations for 
the solar peculiar motion. The classical result, quoted 
in Mihalas and Binney (1981), gives {Uq, Vq, Wq) = 
(+9.0, +12.0, +7.0) km s"-'^, essentially based upon 
Delhaye's (1965) compilation. More recent values 
for the solar peculiar motion do not seem to have 
converged to a single value, especially in the V- 
component. Dehnen and Binney (1998) find {Uq, Vq,Wq) = 
(+10.00 ± 0.36, +5.25 ± 0.62, +7.17 ± 0.38) km s'^ 
from a carefully selected unbiased sample of Hip- 
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parcos stars, while Miyamoto and Zhu (1998) find 
{Uq,Vq,Wq) = (+10.62±0.49,+16.06±1.14,+8.60± 
1.02) km s~^ from 159 Cepheids, also from the Hip- 
parcos catalogue. As suggested earlier by Robin and 
Oblak (1987), such variation in the T^q component is 
probably due to difficulties in separating the asym- 
metric drift from the intrinsic solar motion, and also 
because of the peculiar motions exhibited by the very 
young OB stars, which are still moving under the in- 
fluence of the spiral arm kinematics and/or of their 
parent molecular cloud. 

In 1985, the lAU adopted a value for the motion 
of the LSR of 220 km s"^ Kerr and Lynden-Bell 
(1986) have discussed extensively the determinations 
of V{Ro), as well as Rq, and the Oort constants A 
and B available until then. From a straight mean of 
different determinations they obtained Vi^sRiRo) = 
222 ± 20 km s^^ (their Table 4), while from inde- 
pendent determinations of Rq, A, and B (their Ta- 
bles 3 and 5), we find yLSR(i?o) = 226 ± 44 km s"!. 
Although the uncertainties involved in Vlsr(-Ro) are 
larger than those of Vq , we shall see that the relative 
motion of disk stars is more affected by uncertainties 
in the Solar peculiar motion than uncertainties in the 
motion of the LSR, since the whole nearby disk pop- 
ulation is approximately co-moving with the LSR. 

4. Model comparisons to the absolute proper- 
motions 

The basic standard kinematic parameters employed 
in the model are described in Mendez and van Altena 
(1996) and Mendez et al. (1999). The standard model 
adopts {Uq,Vq,Wq) = (-hll.O,-f 14.0, -1-7.5) km s"! 
for the peculiar solar motion, and a flat rotation curve 
with VhsniRo) = +220 km s~^. Because we are an- 
alyzing data near the SGP, our model predictions do 
not span a large range in Galactocentric distance, and 
therefore the model is insensitive to the value adopted 
for the slope of the rotation curve near Rq (although 
see Mendez et al. 1999). As for the disk kinemat- 
ics, we have adopted the velocity dispersions given in 
Mendez and van Altena (1996), with the scale- heights 
adopted in Mendez et al. (1999) which reproduce the 
observed SPM color- and magnitude counts, a scale- 
length of 3.5 kpc (again, the model is not sensitive 
to this last parameter, as it enters as a function of 
the Galactocentric distance, which for the SPM-SGP 
data is quasi-constant), and a value of g = for the 
velocity ellipsoid (see Eq. 4 in Mendez and van Altena 



1996); q = implies a velocity ellipsoid parallel to the 
Galactic plane at all heights from the Galactic plane, 
while (J = 1 is for a velocity ellipsoid that points to- 
wards the Galactic center). As shown by Mendez et 
al. (1999) the exact choice of kinematical parame- 
ters for the thick-disk and halo are not critical to the 
conclusions of this Letter. The mean velocity of stars 
with different velocity dispersion is computed via the 
proper moments of the CoUisionless Boltzmann equa- 
tion (CBE). The CBE allow us to relate our adopted 
velocity dispersions to the density laws, and to derive 
an expression for the velocity lag as a function of the 
net rotation of the disk (the local LSR) and the pres- 
sure support provided by the velocity dispersions (see, 
e.g., Kuijken and Gilmore 1989, Gilmore et al. 1989). 
We must emphasize however that, if the data contain 
information that is inconsistent with the assumptions 
made (in particular, the validity of an axisymmetric 
CBE, and the neglect of stellar streams), then these 
assumptions inevitably could lead to biased results - 
the classical problem of all parametric methods. 

4.1. The Solar motion and the LSR speed 

Figures (0) and (||) show the model-predicted me- 
dian absolute proper-motions for the standard model 
parameters compared to the observations. These fig- 
ures clearly show that, while the standard model pro- 
duces a very good fit to the U-component of the 
proper- motion (along the Galactocentric direction), 
the V-component (along Galactic rotation) is grossly 
underestimated, especially for Bj < 14, where the 
disk component dominates the overall kinematics. 

A radical change in the model predictions oc- 
curs when we change the peculiar Solar motion from 
the standard value of Vq = -1-14 km s~^ to V© = 
-1-5 km s~^. The largest effect occurs at the bright- 
est bins, i.e., for nearby disk stars where the Solar 
peculiar motion dominates the reflex motion. On the 
other hand, we may keep the conventional Solar pecu- 
liar velocity, but change the overall rotation speed for 
the LSR from the lAU adopted value of -1-220 km s"^ 
to -1-270 km s~^, as suggested by recent Hipparcos re- 
sults (Miyamoto and Zhu 1998). In this case, the 
change in the predicted motion becomes more im- 
portant at fainter magnitudes where one is sampling 
objects from the other Galactic components located 
at larger distances, and where the dominant effect is 
that of the overall rotation of the disk with respect 
to galaxies, and the relative state of rotation between 
the different Galactic components. At the brightest 
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bins, the eflFect of changing Vlsr(-Ro) becomes also 
noticeable because of the larger fraction of bright gi- 
ants, which can be seen to large distances and, hence, 
the differential rotation effects become amplified. A 
model in which we simultaneously change the So- 
lar peculiar motion and the LSR rotational speed to 
-1-5 km s~^ and -1-270 km s~^, respectively, produces 
a much better fit to the observed motions. We note 
that, in this model, the velocity lags for the thick- 
disk and halo are increased in proportion to the in- 
crease of the disk's rotational speed, as the net ro- 
tation of those two components is kept constant at 
180 km s~^ and km s~^, respectively (for details see 
Mendez et al. 1999). We conclude that a model with 
Vq = +5 km and Vlsr(-Ro) = 270 km clearly 
provides a much better fit to the median motion in V 
than does the standard model with Vq = -1-11 km s~^ 
and VisR(^o) = 220 km s^^. 

The predicted median proper-motion in the U- 
component shows a good fit to the observed values, 
and the changes in the V-component described above 
do not affect this parameter substantially because of 
the orthogonality of the projection effects toward the 
Galactic poles. These results do show us, though, 
that the adopted value for the Solar peculiar motion 
in this direction is the correct one. Figure (|^) shows 
the effect of changing the standard Uq = -1-11 km s"^ 
by ±3 km s^^. An eye-ball fit from Figure (|^) sug- 
gests for the U component of the solar motion a value 
of Uq = -1-11.0 ± 1.5 km s^^. Also, there is no in- 
dication of a local expansion or contraction of the 
Galactic disk, as is also found from the kinematics of 
local molecular clouds (Belfort and Crovisier 1984). 

A change of the orientation of the velocity el- 
lipsoids from a cylindrical to a spherical projection 
tends to produce a slightly larger value for the veloc- 
ity lag, especially at fainter magnitudes. The upper 
solid line in Fig. (J^ shows that, while the fit of the 
model to the observed data for Vq = +5 km s""'^, 
Vlsr(^o) — +270 km s~^ is good at bright magni- 
tudes (5,1 < 14), the model underestimates the lag at 
fainter magnitudes. However, by changing the orien- 
tation of the velocity ellipsoid we can actually increase 
the lag. This is shown in Figure ^ where a model 
with Vq = +5 km s"\ ^lsr(^o) = +270 km s~\ and 
q — 1 gives a better overall fit to the observed median 
motion in the V-component than a model with q = 0. 



4.2. Implications for the mass of the Galaxy 

The large LSR speed suggested by the SPM data 
poses an interesting problem to dynamical models 
of the Galaxy. Since the LSR is, by definition, in 
centrifugal equilibrium with the Galactic potential, 
a larger LSR velocity implies a larger mass for the 
Galaxy interior to the Solar circle. As emphasized 
by Rohlfs and Kreitschmann (1988), practically all 
Galactic-mass models follow, to within 10%, the rela- 
tion: 

= 2.3 X 10- f ^U^V (1) 

where M{r) is the mass interior to the radius r, and 
Vrotir) is the circular speed at a distance r from the 
Galactic center. We see that an increase in the LSR 
speed from 220 km s^^ to 270 km s^^ would imply as 
much as 50% more mass interior to the Solar circle. Is 
there any other evidence that this might be plausible? 

An indication that this is perhaps not only possi- 
ble but, indeed, dynamically required, comes from a 
recent measurement of the absolute proper-motions 
of the Large Magellanic Clouds (LMC) by Anguita 
et al. (1999). They have used QSOs as an inertial 
reference frame to obtain absolute proper-motions of 
stars in three fields toward the LMC. The proper- 
motions are based on CCD frames taken at 7 to 
12 epochs spanning eight years, with a mean accu- 
racy for the proper-motions from these three fields 
of ±0.2 mas yr~^. Their derived motion for the 
LMC is larger than that inferred by previous pho- 
tographic studies that have used either bright stars 
or extended galaxies as a reference frame. Anguita et 
al. derive a mass for the Milky Way out to 50 kpc of 
3 — 4x 1O"'^^M0, in contrast to previous determinations 
(e.g., Lin and Lynden-Bell 1982) that implied a value 
of 7 X 10^^ M0. Anguita's result is fully consistent 
with our findings: Galactic mass models predict that 
the mass interior to the Solar circle is only about 4% 
of the total Galactic mass (Rohlfs and Kreitschmann 
1988). The range of mass values found by Anguita 
et al. implies, through Eq. (1) for r = 8.5 kpc, a 
rotational velocity at the Solar Galactocentric dis- 
tance in the range 250 — 286 km s^^, which indeed 
brackets our results. Interestingly enough, a mass 
of 7 X 10'^"'^ M0 implies a very improbable rotational 
speed of 120 km s~^. 

On the other hand, our derived value for Vlsr(^o) 
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seems to be in conflict with a recent determination of 
tfie proper motion of Sgr A* by Reid et al. (1999), and 
Backer & Sramek (1999) who find VLSR(i?o)/i?o « 
29 kms^^kpc^^ (employing the good assumption 
that Sgr A* is at rest), which would require Rq larger 
than 9 kpc, inconsistent with almost all previous de- 
terminations. Could this inconsistency point, per- 
haps, to evidence against our basic model assump- 
tions which would then lead to a biased I^lsr(^o)? 

4.3. Other evidence for a larger Galactic mass 

Recently, Zaritsky (1999) has reviewed most of the 
evidence regarding the mass of the Galactic halo and, 
by extension, the total mass of the Galaxy. He points 
out that a squared ratio of rotational circular ve- 
locities of 1.5 (such as that found from the 220 to 
270 km s~^) is not considered a serious discrepancy 
in the modeling of the total Galactic mass. He fur- 
thermore points out that timing arguments involving 
the relative motion of M 31 and the Milky Way imply 
a mass for the latter of ~ 1.4 x 10^^ M©. A simi- 
lar argument, based on the Leo I-Milky-Way system, 
leads to a mass of 1.1 — 1.5 x 10"'^^ Mq. In both cases, 
the analysis excludes angular momentum, overlapping 
mass distributions at earlier times, and the evolution 
(growth) of the Galactic halo with age, all of which 
would lead to an increased total mass. Zaritsky also 
reviews an extension of the two-body calculations into 
the larger environment, which leads to an even larger 
mass, in the range 1.9-2.3x10^2 Mq, concluding that 
all this evidence implies a mass out to ~ 200 kpc of 
> 1.2 X 10^2 Mq. Thus, Anguita's results based on 
the LMC proper motions, and its correspondence to 
our findings, are in line with Zaritsky's arguments. 

5. Conclusions 

We have used the distribution of absolute proper 
motions of ^ 30, 000 randomly selected stars, part 
of the Southern Proper Motion survey, to constrain 
a structural & kinematic model of the Galaxy. The 
absolute proper-motions in the U-componcnt indicate 
a solar peculiar motion of 11.0 ± 1.5 km s^^, with no 
need for a local expansion or contraction term. In the 
V-component, the absolute proper- motions can only 
be satisfactorily reproduced by the model if we adopt 
a solar peculiar motion of -t-5 km s~^, a large LSR 
speed of 270 km s~^, and a (disk) velocity ellipsoid 
that always points towards the Galactic center. 

The larger than expected LSR speed leads to a 



larger mass for the Galaxy interior to the Solar Circle 
by a factor of about 1.5. This is compared with other 
recent evidence indicating that the Milky- Way might 
indeed be more massive than previously thought. 
This evidence comes mainly from a new measure- 
ment of the proper motion of the LMC (Anguita et 
al. 1999) and from the binary motion ("timing argu- 
ment" ) of the M31 - Milky Way and Leo I - Milky Way 
pairs (Zaritsky 1999). Our larger LSR speed is also 
coincident with a similar finding by Miyamoto and 
Zhu (1998) using 0-B5 stars and Cepheids from the 
Hipparcos catalogue. We stress that our result of a 
large LSR speed rests crucially on the assumptions of 
both axisymmetry and equilibrium dynamics. Failure 
to meet these conditions would lead to inconsisten- 
cies, such as those related to the Solar-galactocentric 
distance indicated at the end of Section 4.2. 

High accuracy proper motion measurements for the 
SMC, and proper motions for some of the nearby 
dwarf Spheroidals will help in confirming a larger 
value for the mass of the Galaxy, while future space 



interferometric space missions like FAME (http : //www . usno . navy 



GAIA (^tp : //astro . estec . esa . nl/GAIA/| ), and SIM 
( b.ttp : //sim. jpl .nasa. gov/ ), will directly measure 
the disk rotational speed across the entire Galaxy. 
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Fig. 1 . Median observed (fillc;d dots with error bars) 
and predicted (lines) absolute proper-motions along 
U (Galactocentric direction). In some cases the er- 
ror bars are smaller than the size of the dot. The 
solid line is for the standard model, the dashed line 
is for a model with a Solar peculiar motion in the U- 
componcnt of +8 km s^^, while the dot-dashed line 
is for a model with Uq = +14 km s~^. The stan- 
dard model, with Uq = +11 km s~^ provides the 
best fit to the U-componcnt of the observed motion. 
The agreement to the model also indicates the ab- 
sence of any significant expansion/contraction in the 
local disk (see text). 

Fig. 2. — Median observed (filled dots with error 
bars) and predicted (lines) absolute proper-motions 
along V (Galactic rotation). In some cases the er- 
ror bars are smaller than the size of the dot. The 
lower solid line is from the standard model with 
Vq = +U km s-\ and yLSR(i?o) = +220 km s'^ 
The dashed line is for a model with a disk having a 
rotational speed of -1-270 km s~^. The dotted line 
indicates the predictions for a model with a Solar 
peculiar motion in the V-component of -1-5 km s~^ 
instead of the classical value +14 km s~^ adopted 
in the standard model, while the triple-dot dash line 
shows a model with, both, Vq = +5 km s"-'^, and 
VLSR(i?o) = +270 km s'^ Clearly, a larger LSR 
speed improves the model predictions at the fainter 
{Bj > 15) bins. Finally, the upper solid line is 
for a model with Vq = +5 km s^^, VisR(-Ro) = 
+270 km s~^, and q = 1, indicating that the best fit 
to the observed median motion in the V-component 
is provided by a velocity ellipsoid whose major axis 
points to the Galactic center as we move away from 
the Galactic plane. 
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